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Abstract 



Experimental data indicate small spin-orbit splittings in hadrons. 
For heavy-light mesons we identify a relativistic symmetry that sup- 
presses these splittings. We suggest an experimental test in electron- 
positron annihilation. Furthermore, we argue that the dynamics neces- 
sary for this symmetry are possible in QCD. 
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I. INTRODUCTION 



Recently, Isgur has re-emphasized the experimental fact that spin-orbit split- 
tings in meson and baryon systems, which might be expected to originate from one- 
gluon-exchange (OGE) effects between quarks, are absent from the observed spec- 
trum. He conjectures that this is due to a fairly precise, but accidental, cancellation 
between OGE and Thomas precession effects, each of which has "splittings of hun- 
dreds of MeV" I]. Taking the point of view that precise cancellations reflect symme- 
tries rather than accidents, we have examined what dynamical requirements would 
lead to such a result. One of us recently observed |§ that a relativistic symmetry is 
the origin of pseudospin degeneracies first observed in nuclei more than thirty years 
ago [§,|J. We find that a close relative of that dynamics can account for the spin 
degeneracies observed in hadrons composed of one light quark (antiquark) and one 
heavy antiquark (quark). 

Below, we first elucidate the experimental evidence for small spin-orbit splittings. 
Then we identify the symmetry involved in terms of potentials in the Dirac Hamil- 
tonian for heavy-light quark systems, and note the relation to the symmetry for 
pseudospin. We show that the former symmetry predicts that the Dirac momentum 
space wavefunctions will be identical for the two states in the doublet, leading to a 
proposed experimental test. Finally, we argue that the required relation between the 
potentials may be plausible from known features of QCD. 

II. EXPERIMENTAL AND LATTICE QCD SPECTRUM 

In the limit where the heavy (anti) quark is infinitely heavy, the angular momen- 
tum of the light degrees of freedom, j, is separately conserved |J. The states can 
be labelled by lj, where I is the orbital angular momentum of the light degrees of 
freedom. In non-relativistic models of conventional mesons the splitting between l l+ i 
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and levels, e.g. the p3 and pi or ds and dz levels, can only arise from spin-orbit 
interactions The pi level corresponds to two degenerate broad states with differ- 
ent total angular momenta J = j ±sq (here j = |), where sq is the spin of the heavy 
(anti)quark ||. For example, in the case of Z)-mesons, Sq = | and the two states 
are called Dq and D[. There are also two degenerate narrow p-i states D\ and D* 2 
||. The degenerate states separate as one moves slightly away from the heavy quark 
limit, and their spin- averaged mass remains approximately equal to the mass before 
separation. 

For the D-mesons, the CLEO collaboration claims a broad J p = 1 + state at 
2461±34 ± 10 ± 32 MeV @, belonging to the pi level, in close vicinity to the D* 2 
at 2459 ± 2 MeV ]7[|, belonging to the ps level, indicating a remarkable ps-pi spin- 
orbit degeneracy of —2 ± 50 MeV. It is appropriate to extract the spin-orbit splitting 
this way since: Firstly, the charm quark behaves like a heavy quark. Secondly, 
the difference between the D[ and D\ levels is the best indicator || of the pa-pi 
splitting in the absence of experimental dataQ on the D$, as opposed to the difference 
between the D[ and spin-averaged p-i level at 2446 ± 2 MeV. Spin-averaged masses 
are determined from experiment . 

For the iT-mesons, the pi level is at 1409 ± 5 MeV, with pa nearby at 1371 ± 3 
MeV, corresponding to a pz-pi splitting of —38 ± 6 MeV. The splitting between the 
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higher-lying ds and ds levels is —4 ± 14 MeV or 41 ± 13 MeV, depending on how the 
states are paired into doublets. These results indicate a near spin-orbit degeneracy 
if the strange quark can be treated as heavy, although it has certainly not been 
established that such a treatment is valid. 

For 5-mesons, both L3 |5| and OPAL M have performed analyses, using input 



*The FOCUS collaboration preliminarily found Dq at a mass of 2420 MeV [31|. The error 
on the mass was not reported. 
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from theoretical models and heavy quark effective theory, to determine that the ps-pi 
splitting is 97 ± 11 MeV (L3) or -109 ± 14 MeV (OPAL). Note that these are not 
model-independent experimental results. In the same analyses the mass difference 
between the B% and Bq, an approximate indicator of the ps-pi splitting, is 110 ± 11 
MeV (L3) or -89 ±14 MeV (OPAL). The L3 result agrees with lattice QCD estimates 



of 155i? 3 ±32 MeV |10| and 183±34 MeV fil l. However, according to other estimates 



T2| , the splitting is less than 100 MeV, and consistent with zero. Recently, 31 ± 18 



MeV was calculated [T3|]. One lattice QCD study found evidence for a change of 
sign in the splitting somewhere between the charm and bottom quark masses, albeit 



with large error bars fll4| . A splitting of 40 MeV serves as a typical example of 
model predictions ||15|| , although there is variation in the range -155 to 72 MeV [ |I~6|j , 
summarized in ref. JT3[ . 

In order to more quantitatively measure the spin-orbit splitting, define 

(pa - pi) 

r = — * (1) 

((4 P| ±2 P| )/6- S i)' 1 J 

where all entries refer to masses. The experimental data on D, K and B mesons 
give respectively r = 0.00 ± 0.10, -0.06 ± 0.00 and 0.23 ± 0.04 (L3) or -0.23 ± 0.03 
(OPAL). For the Dirac equation with arbitrary vector and scalar Coulomb potentials, 
the only cases for which the relevant analytic solutions are known, —0.7 ~ r ~ 0.6. It 
is hence evident that the spin-orbit splittings extracted from experimental results are 
indeed small. 

There is also evidence in light quark mesons and baryonic systems that the spin- 
orbit interaction is small [I]]. In non-relativistic models, meson and "two-body" 
baryon spin-orbit interactions are related and, for a specific class of baryons, the 
spin-orbit interaction is small for exactly the same reasons that it is small in mesons 

iifl. 
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III. A DYNAMICAL SYMMETRY FOR THE DIRAC HAMILTONIAN 



If we consider a system of a (sufficiently) heavy antiquark (quark) and light quark 
(antiquark), the dynamics may well be represented by the motion of the light quark 
(antiquark) in a fixed potential provided by the heavy antiquark (quark). Let us as- 
sume that both vector and scalar potentials are present. Then the Dirac Hamiltonian 
describing the motion of the light quark is 

H = a ■ p + (3(m + V s ) + V v + M, (2) 

where we have set % = c = 1, a, /3 are the usual Dirac matrices, p is the three- 
momentum, m is the mass of the light quark and M is the mass of the heavy quark. 

This one quark Dirac Hamiltonian follows from the two-body Bethe-Salpeter equa- 
tion in the equal time approximation, the spectator (Gross) equation with a simple 
kernel, and a two quark Dirac equation, in the limit that M is large |Jl7Hl9j. If the 
vector potential, Vv(r), is equal to the scalar potential plus a constant potential, U, 
which is independent of the spatial location of the light quark relative to the heavy 
one, i.e., Vy(f) = Vs(r) + U, then the Dirac Hamiltonian is invariant under a spin 
symmetry [pO| , pT |, [H , Si] = 0, where the generators of that symmetry are given by, 
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where Sj = fXj/2 are the usual spin generators, Uj the Pauli matrices, and Sj = U p §i U p 
with U p = —. Thus Dirac eigenstates can be labeled by the orientation of the spin, 
even though the system may be highly relativistic, and the eigenstates with different 
spin orientation will be degenerate. 

For spherically symmetric potentials, Vy(f) = Vy(r), Vs(r) — V s (r), the Dirac 
Hamiltonian has an additional invariant algebra; namely, the orbital angular momen- 
tum, 



where t{ = U p £i U p and £i = (fx This means that the Dirac eigenstates can 
be labeled with orbital angular momentum as well as spin, and the states with the 
same orbital angular momentum projection will be degenerate. Thus, for example, 
the n r pi/2 and n r p 3 / 2 states will be degenerate, where n r is the radial quantum 
number. 

Thus, we have identified a symmetry in the heavy-light quark system which pro- 
duces spin-orbit degeneracies independent of the details of the potential. If this 
potential is strong, the heavy-light quark system will be very relativistic; that is, the 
lower component for the light quark will be comparable in magnitude to the upper 
component of the light quark. It is remarkable that non-relativistic behaviour of 
energy levels can arise for such fully relativistic systems. 

This symmetry is similar to the relativistic symmetry identified as being re- 
sponsible for pseudospin degeneracies observed in nuclei [§,|J. In contrast to spin 
symmetry, pseudospin symmetry has the pairs of states ((n r — l)si/ 2 , ^^3/2), 
((n r — l)p 3 /2, n r /5/2), etc. degenerate, making the origin of this symmetry less trans- 
parent. The pseudospin generators are 
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For pseudospin symmetry, the nuclear mean scalar and vector potential must be equal 
in magnitude and opposite in sign, up to a constant, Vy = —Vs + U . Relativistic 
mean field representations of the nuclear potential do have this property; that is, 



Vs ~ — Vy [P2| , p3| . We will return later to the question of whether the relation 
V v = V s + U arises in QCD. 

It has previously been observed that pseudospin symmetry improves with in- 
creasing energy of the states, for various potentials 0. A similar behaviour may 
be expected for spin symmetry, consistent with the experimental observations that 
spin-orbit splittings decrease for higher mass states 
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The Dirac Hamiltonian (fj) encompasses the effects of the OGE and Thomas pre- 
cession spin-dependent terms customarily included in non-relativistic models ]![. 

IV. EXPERIMENTAL TEST 

In the spin symmetry limit, the radial wavefunctions of the upper components of 
the Dirac wavefunction of the two states in the spin doublet will be identical, be- 
having "non-relativistically" , whereas the lower components will have different radial 
wavefunctions. This follows from the form of the spin generators given in Equation 
(0). The (1, 1) entry of the operator matrix is simply the non-relativistic spin opera- 
tor which relates the upper component of the Dirac wavefunction of one state in the 
doublet to the upper component of the other state in the doublet. Since this operator 
does not affect the radial wavefunction, the two radial wavefunctions must be the 
same. By contrast, the lower component wavefunction is operated on by U p which 
does operate on the radial wavefunction because of the momentum operator. 

As an example, we show in Figure 1 the upper and lower components for Dirac 
wavefunctions of the pi/2 — P3/2 doublet. The scalar and vector potentials were deter- 
mined by matching the available spectral data of the D-mesons, assuming a ps — pi 
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splitting at the lower end of the range defined by the experimental value of —2 ± 50 
MeV. This maximizes the wavefunction differences. In this realistic case, Vy ~ Vs+U, 
so the radial wavefunctions for the upper components are not exactly identical but 
are very close, whereas the radial wavefunctions for the lower components are very 
different. 

Likewise the momentum space wavefunctions for the upper components will be 
very similar, as seen in Figure 2, again because the spin operator does not affect the 
wavefunction. However, since U p depends only on the angular part of the momentum, 
p = 2 , it does not affect the radial momentum space wavefunction. In Figure 2 we 
see that the radial momentum space wavefunctions are very similar for the lower 
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components as well. This prediction of the symmetry can be tested in the following 
experiment. 

The annihilation e + e" — > D^Dq, D^D 2 and D 2 D 2 allows for the extraction of 
the Dq and D 2 electromagnetic static form factors and the Dq to D 2 electromagnetic 
transition form factor. The photon interaction ensures that all radial wavefunctions 
of the light quark are accessed. When spin symmetry is realised, there are only 
two independent radial momentum space wavefunctions, which should enable the 
prediction of one of the three form factors in terms of the other two. This should 
enable the verification of the predictions of spin symmetry. On the other hand, non- 
relativistic models, with no lower components for the wavefunctions, have only one 
independent radial wavefunction, which will lead to the prediction of two of the form 
factors in terms of the remaining one. This might be too restrictive. The proposed 
experiment can be carried out at the Beijing Electron Positron Collider at an energy 
of approximately 1 GeV above the ^(4040) peak in the final state DDirir. 

An equivalent experiment for K-mesons would involve detection of the KKtttt 



final state, which has already been measured p4|. The wavefunctions of i^-mesons 
fitting the experimental spectrum show similar behaviour to the D-mesons, with the 
C3 and pi wavefunctions even more similar than in Figures la and 2. 

2 2 

If B-mesons do also exhibit spin symmetry, one can do equivalent experiments 
around 1 GeV above the T(3S) peak at the SLAC, KEK or CESR B-factories. 

V. QCD ORIGINS 

If such a dynamical symmetry can explain the suppression of spin-orbit splitting in 
the hadron spectrum, the question remains as to why it might be expected to appear in 
QCD. To address this, we first recall the ongoing argument as to whether confinement 
corresponds to a vector or scalar potential |25 |. The first natural expectation was that 



confinement reflected the infrared growth of the QCD coupling constant, enhancing 



the color-Coulomb interaction at large distances, see e.g. Ref.( ||26|| ). An involved 
two- (or multi-) gluon effect has been proposed |27j to account for the origin of a 
scalar confining potential. 

The existence of one or the other of these vector and scalar potentials is not neces- 
sarily exclusionary - they may both be realised. The arguments in Ref.( suggest 
further that they are related, with the scalar exceeding the vector by an amount which 
may be approximately constant as one saturates into the linear confining region at 
large separations. We very briefly reiterate the basic argument of Ref.( |28fl ) here. 

The starting point is to accept the standard approach |26] that renormalization- 
group-improved single-gluon-exchange produces a linearly increasing vector potential 
between a quark and an antiquark. One then considers what to expect for multiple 
gluon exchange, starting with two gluons. Since two gluons are attracted to each other 
in a color singlet channel, and also have a zero mass threshold (as for massless quark- 
antiquark pairs), it is reasonable to conclude that a (Lorentz and color) scalar gluonic 
condensate develops, along with a mass gap for a glueball state. These developments 
are indeed observed in lattice QCD calculations. 



Ref.( p8|) goes on to argue that renormalization-group-improved single-glueball- 
exchange involves the square of the QCD coupling and so, despite the massiveness 
of the object exchanged, also leads to a (now scalar) confining potential between 
quarks and antiquarks. This further implies that the ratio of the slopes of the two 
potentials in their common linear (confining) region is given by the square of the ratio 
of the QCD scale for growth of the coupling constant to the value of the mass gap 
of the condensate formation. This ratio may be expected to be of order one as both 
quantities are determined by the underlying QCD scale. 

If the two potentials do indeed have similar slopes in the region outside that 
dominated by the color Coulomb interaction, they would necessarily differ only by 
an approximately constant value, in that region. Thus, the origin of the dynamical 
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symmetry may not be unreasonable, and may indeed be a natural outcome of non- 
perturbative QCD. 

On the other hand, identically equal vector and scalar potentials, except for a 
constant difference, would appear to be coincidental. An ameliorating effect is that 
to produce an approximation to the spin symmetry of Eq. (2) this condition need 
only hold in regions where the wavefunctions are substantial. 

The determination of QCD potentials, from models like the minimal area law, 
stochastic vacuum model, or dual QCD, and from lattice QCD, is hampered by the 
problem of rigorously defining the concept of a potential from QCD when one quark 
is light. It suffices to say that there is no agreement on the mixed Lorentz character 



of the potential even between two heavy quarks ||29|| , where the potential can be 
rigorously defined, although lattice QCD results are consistent with simply a vector 
Coulomb and scalar linear potential [RO . 



VI. SUMMARY 

The observation of "accidental" spin-orbit degeneracies observed in heavy-light 
quark mesons can be explained by a relativistic symmetry of the Dirac Hamiltonian 
which occurs when the vector and scalar potentials exerted on the light quark by 
the heavy antiquark differ approximately by a constant, Vy ~ V$ + U. Conversely, 
if future experiments determine that spin-orbit splittings are small not only for the 
lowest excited states in mesons but are small throughout the meson spectrum, this 
experimental fact dictates that the effective QCD vector and scalar potentials between 
a quark and antiquark are approximately equal up to a constant, which would be 
a significant observation about the nature of non-perturbative QCD. Furthermore, 
the approximate symmetry predicts that the spatial Dirac wavefunction for the spin 
doublets will be approximately equal in momentum space, a feature which can be 
tested in electron-positron annihilation. We have argued that Vy ~ Vg + U may occur 
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in QCD, particularly for regions of space dominated by the light quark wavefunction. 
Work is in progress to extend this symmetry to purely light quark systems. 

This research is supported by the Department of Energy under contract W-7405- 
ENG-36. 
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FIGURES 



Figure 1: (a) The square of the Dirac radial wavefunction of the upper component 
times r 2 . (b) The square of the Dirac radial wavefunction of the lower component 
times r 2 . pa is the solid line and m is the dashed line. Note that the lower component 

2 2 

is comparable to the upper component. The wavefunctions are solutions of the Dirac 
equation (see Eq. (§)) with Coulomb potentials Vs(r) = ^f+C/s and Vv(r) = ^-+Uv, 
where a s = -1.279, U s = 506 MeV, a v = -0.779, U v = 515 MeV, m = 330 MeV 
and M = 1480 MeV. This corresponds to a pa — pi splitting of -52 MeV. 

Figure 2: (a) The square of the Dirac momentum space wavefunction of the upper 
component times q 2 . (b) The square of the Dirac momentum space wavefunction of 
the lower component times q 2 . Other conventions are the same as in Figure 1. 
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